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We generalize Hedin’s equations to a system of superconducting electrons coupled with a system 
of phonons. The electrons are described by an electronic Pauli Hamiltonian which includes the 
Coulomb interaction among electrons and an external vector and scalar potential. We derive the 
continuity equation in the presence of the superconducting condensate and point out how to cast 
vertex corrections in the form of a non-local effective interaction that can be used to describe both 
fluctuations of spin and superconducting phase beyond the screened Coulomb self-energy diagram. 


I. INTRODUCTION 

Hedin’s formally exact iterative procedure to generate 
the one particle Green’s function for electrons that in¬ 
teract via the Coulomb potential 1 has provided a useful 
basis for further approximations. Especially the famous 
GW approximation, as already introduced in the original 
paper, is very successful in the first principle analysis of 
physical system using Many-Body-Perturbation theory 2 . 
Nowadays, with increasing computational power, even a 
fully self-consistent treatment of the screened Coulomb 
potential W within Hedin’s cycle for both, solids and 
molecules is in reach 3,4 . 

While we thus have a firm foundation of ab-initio 
electronic structure perturbation theory, to calculate 
superconductivity (SC), the Green’s function approach 
boils down to the Eliashberg equations 5- '. There, the 
Coulomb potential is usually reduced into the /x* pseudo- 
potential 6, ' and the phonons are computed externally, 
i.e. in the normal state. This is reasonable since the vi¬ 
brational structure and the renormalized Coulomb po¬ 
tential appear to be largely independent on the SC con¬ 
densation. Thus, in the usual Eliashberg approach, SC 
is computed starting from a converged normal state elec¬ 
tronic band and phonon structure 6 . On the other hand, 
spin-fluctuations are discussed as a pairing mechanism 
in the unconventional SC such as the cuprates and iron 
pnictides 8 . With fluctuations we mean here self-energy 
corrections that describe effects beyond the screened 
Coulomb diagram. The paramagnon peak as the quasi¬ 
particle that corresponds to the spin-fluctuation inter¬ 
action becomes entirely different in the SC phase 9 and 
it therefore is essential to consider SC and the screen¬ 
ing of electronic interactions on the same footing. The 
paramagnetic spin-fluctuations are included in ab-initio 
methods on a semi-ab-initio level so far. Essenberger et 
al. 10 have derived an effective interaction that is based 
on a self-energy contribution due to the magnetic sus¬ 
ceptibility and the exchange-correlation kernel of time 
dependent density functional theory 11 . We note, how¬ 
ever, that the effective interaction derived in Ref. 10 is in 
the normal state and the self-energy that is used may 
suffer from double counting. The diagrams considered 
for the construction of the effective interaction are for¬ 
mally similar to approach of Berk and Schrieffer 12 . In 
the context of a Hubbard model, this approach has been 


extensivley used to describe the SC phase of the iron 
based superconductors (see Ref. 13 for a review). 

In a similar manner fluctuations of the SC order 
parameter 14 in constrained geometries are not well de¬ 
scribed in present ab-initio methods. While the famous 
Mermin-Wagner theorem 15 forbids the ordering in 2D 
due to the onset of long wavelength fluctuations it is still 
unclear to what extend this result has implications in ex¬ 
periment. For example, SC has been observed in a lead 
surface on a silicone substrate down to the single atomic 
layer 16 . Similarly, taking the Coulomb interaction into 
account, the excitation energies of the collective states 
of paired electrons starts at the plasma frequency 17-19 . 
Since plasmons with a small excitation energy may ex¬ 
ist in low dimensional systems, the collective excitations, 
i.e. SC fluctuations, may have to be taken into account 
self-cons istently. 

To present a fundamental theory to describe SC based 
on the basic interactions among electrons is the content 
of this paper. In particular, we derive Hedin’s equations 
for SC electrons that interact via the Coulomb poten¬ 
tial and with a system of phonons in Sec. III. Our ap¬ 
proach treats the electronic screening, the coupling to 
the phonons and SC on the same footing. In Sec. V, we 
point how to interpret vertex corrections as an effective 
interaction, without introducing double counting with re¬ 
spect to the screened Coulomb interaction. Our coupled 
equations are capable to describe fluctuation effects in 
principle exactly. This means in particular spin angular 
momentum transfer processes in a paramagnetic system 
as well as, in the language of Anderson 1 ', iso-spin mo¬ 
mentum transfer with respect to the Nambu off diagonal 
components, which corresponds to the Cooper-channel, 
of the Green’s function. 


II. THE HAMILTONIAN 

In this Section we define the basic Hamiltonian govern¬ 
ing our system of electrons and phonons. Before we intro¬ 
duce the individual parts of the Hamiltonian in Eq. (5), 
we define the notation that we use in this work in the 
following separate paragraph. 

a. Notation To simplify the notation we 
introduce the usual thermal average (•••)o = 

Tr{e _ ^ J7_/iJ ^...}/Tr{e _ ^ 77_/iJ ^} and (A) T = 


2 


(TU((3,0)A) 0 /(U(f3,0))o with the time ordering symbol 
T and the time evolution operator U(t, t'). 

Following the notation of Nambu 18 and Anderson 17 , in 
order to describe it as a usual single particle Hamiltonian 
we introduce the Nambu field via 

^( r ) = ( i>{r t) i>(r J.) V’ t ( r t) . (1) 

We use T ns to transpose in Nambu (n) and spin space (s) 
and label Nambu components with a = 1,-1 and spin 
components with p =f, J,. Tr n3 means to take the trace 
in Nambu and spin space. 

We further group variables in the notation 

1 = rin 1 = ai/xi 

fdl = /dn/drr = E wai / 2 ) 

§i ,2 = 6(ri-r 2 )x bj 2 = 6 Ql Q 2 6 Ml/i2 
x6(n - r 2 ) 

The bar indicates Nambu-spin. Thus, similarly, we in¬ 
dicate matrices in Nambu and spin space with a bar on 
top. For example, in this notation, the electronic single 
particle Green’s function reads 

G(rr, r'r') = -(^(rr) ® ( 7 'V)) T . (3) 

The symbol ® means to take the outer product in Nambu 
and spin space and the subscript I refers to an opera¬ 
tor in the Heisenberg picture 0\(t) = U(0,t)OU(t,0). 
We promote the Nambu and spin indices to the argu¬ 
ment to refer to the matrix elements, similar to the 
usual matrix notation. There, for a matrix A, one refers 
to the matrix elements as A 7J . 0(11,22) is thus the 
aiO’i,a 2 o ’2 matrix element in Nambu and spin space - 
a function of two space and time variables 1 and 2, re¬ 
spectively. We use this matrix element notation also for 
higher order tensors, such as 1(11 , 22, 33, 44) which is 
the (ai< 7 i, a 2 02 , 0 : 303 , 0404 ) component of the 4x4 x 4x4 
object 1(1,2, 3,4)- 

If indexes appear on one side of an equation and not on 
the other we define that these are summed or integrated. 

We use the Pauli matrices both in spin as well as in 
Nambu space. For the spin-Pauli matrices we use <To,x,y, z 
while for the latter we use 7 o lX ,y, 2 - 

We introduce the Nambu time ordering symbol 

Ti^i (rr) ® <Pj (r'r') = 0 (t — r')<^i(rr) ® )Pj( r ' T ') 

— 9(t' — t) (^/(r'r 7 ) ® ^(rr)) 1 ' 311 

(4) 

This form is equivalent to a time ordering in every com¬ 
ponent of the Nambu and spin matrix. Note a peculiarity 
in the equal time limit: We define that the a, a 1 = — 1, — 1 
component behaves different to the 1,1 component in the 
sense that here the second argument is taken infinitesi¬ 
mally after the first. The reason is that in the equal 
time limit it is necessary to recover the density opera¬ 
tor n(r) = i/4(r) • ctq ■ V’( r )> also in the Nambu —1, —1 
channel. This equal time limit appears in the interac¬ 
tions in a diagrammatic expansion so the definition must 


be chosen to recover this limit 20 . The missing reorder¬ 
ing of the operators causes an additional minus sign and 
thus ni(rr) = — (l/2)Tr ns {lim T /_>. T r ^ •T!^ ^ I(rr)®!F/(rr , )} 
with t z instead of tq without the reordering. 

b. Contributions to the Hamiltonian Our starting 
point is a system of interacting electrons that is coupled 
with a system of non-interacting phonons 

H = H 0 + H P + H B _ e + H e _ p + ff aux . (5) 

Here we distinguish the single particle part H 0 , the non¬ 
interacting phonon part H p , the electron-electron inter¬ 
action H e _ B , the electron-phonon interaction H 0 _ p and an 
auxiliary external potential that is set to zero after the 
derivation. In the single particle part H 0 we distinguish 
H 0 = H„ + H with a normal state //,, and SC part II . 
The normal state part is H n = f drift'(r) ■ H n (r ) • ift(r) 
and ift(rp) is the electron spin p field operator. Then, 
the normal state contribution of the single particle part 
of the Hamiltonian reads 

H n (r) = er 0 (^(—iV + A ext (r )) 2 + 0 „ t (r)) + 

+ S • £? ext (r). (6) 

Here B oxt (r) = V x A ext (r), S = ( cr x o y a z ) T and 
A oxt (r) is an external vector potential. Further 0 ext (r) 
is the scalar potential including the Coulomb attrac¬ 
tion of the ions in their equilibrium position. The ex¬ 
istence of a pair condensate with a macroscopic num¬ 
ber of electron pairs is measured with the order pa¬ 
rameter of SC x( r ■> r ') = (V , ( r ) • 3? • ift( r ')) where 
ift(r) = ( ift(r f) ift(r j.) ) . The 4 component vector 
$ = ( iciy — a z do &x ) T parametrizes the order param¬ 
eter into 1 singlet dq and 3 triplet parts <4> 2 ,4*31 4 > 4 . The 
singlet and triplet SC part of the total Hamiltonian that 
couples to an external singlet and triplet pairing field 

A LA r , r ') is 



- ft(r) ■ $ • ^(r') ■ A^ t (r, r')) . (7) 

Here, the combination A Bxt (r, r') ■ = -A KXt (r', r) ■ 

is totally antisymmetric. Thus, the pair Hamilto¬ 
nian Eq. 7 is hermitian by construction. We use q = q, X 
with the Bloch vector q and the mode number A to in¬ 
dicate the quantum number on the phonon operators b q . 
Then, the phononic single particle Hamiltonian is 

H p = Y, Q MK + \)- ( 8 ) 

9 

We assume the electrons to interact via a spontaneous 
Coulomb interaction 



x ift\rp)ift\r'p')ift(r'p')ift(rp), (10) 
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with w 0 (r,r') = 1 /|r — r'\. w 0 (r,r') is independent 

of spin and also the screened potential will share this 
property. Following Ref . 21 one could start with a spin 
dependent interaction. In this work we take a different 
route. After the derivation of Hedin’s equations for a SC, 
we discuss the important self-energy diagrams that form 
an effective spin dependent interaction.and the electron 
phonon interaction is (we use the notation —q = — q,X 
and i = 0 , x , y , z) 

H„- p = Y f dr ft( r ) ' ' ^(r)g z (rq)(b q + bl q ) .(11) 

iq 


Because H P _ P is hermitian, gi(rq) satisfies gi(rq) = 
g*(r,—q). The operator H B is not diagonal in the elec¬ 
tronic field operator ip(r). To derive the Hedin equations 
we use the auxiliary external fields q)(rr) and J q (r) that 
are set to zero after the derivation. We define 

H am = H° +H P , (12) 

and choose H® ux (t i) = J driip(l)ijj^ (r 3 ) ■ a 3 • ^(tt) and 
H P U1C = J q (b q +tf_ q ) where J _ 9 = J* to make the operator 
hermitian. 

III. HEDIN’S EQUATIONS FOR 
SUPERCONDUCTORS 

In this Section, we derive the Hedin equations for a 
superconductor. We point out that we can work in the 
ground state, the equilibrium finite temperature 20 or the 
more general Keldysh formalism 22 . The Hedin equation 
formally have the same shape. Here we work in the 
equilibrium finite temperature formalism for definiteness. 
Our starting point is the equation of motion of the sin¬ 
gle particle Green’s function. We express the two par¬ 
ticle part as a functional derivative with respect to the 
auxiliary field. Using a similar strategy for the appear¬ 
ing self-energy and vertex terms leads to a set of closed 
equations, independent on the auxiliary potential - the 
Hedin equations. 


A. Green’s Function Dyson Equation 

In the following we derive the Dyson equation for the 
electronic Green’s function from the equation of motion 

— d Tl G(l 1 , 22) = bq, 2 ^ 1,2 A 

+ ([H,*(r 1 l)] l (n)*l(22)) T .( 13) 

Similarly to the original derivation, we are describing the 
two particle part of the equation of motion as the func¬ 
tional derivative with respect to the auxiliary field <p(l). 
Note here that although H and bH/8ip(l) do not com¬ 
mute in the case of a SC, the time ordering still allows 
us to use j^TU(^,0) = \TU(/3,0)&}(1) • a 0 r z •&i(l). 


We complete the derivation of the Dyson equation in the 
Appendix A with the result 

G(l,2) = G a (l,2) + G H (l,3)-£(3,4)-G(4,2). 

(14) 

Also the self-energy £(3,4) is defined in the Appendix 
A. Our Dyson equation starts from the Hartree Green’s 
function G H (3, 2) that satisfies 

(-T 0 <j 0 d Tl 5!,3 - H h (1 , 3)) ■ G h (3, 2) = Sr, 2 0 - o r 0 

(15) 

where the Nambu and spin matrix opera¬ 
tor H m (l , 3) in first quantization is given by 
H h ( 1,3) = H 0 (l,3) + <P(l)a 0 T z 5 lt3 with the Hartree 
field <P(1) = J n(r)w 0 (r i, r)dr + ip(l) where n(r) is the 
electron density. The bare Hamiltonian is 


H 0 (l,3) 


H»(l)6i, s A^(l,3)-* \ 

-A:J1,3)-* ~(H a (l)) Ts Sr, 3 J ’ 

(16) 


with the normal state part H a of the total Hamiltonian 
given in Eq. ( 6 ). Note that because of the total anti¬ 
symmetry of A kxI (1 , 3) ■ and S Ta = S* the bare Hamil¬ 
tonian H 0 is hermitian. The self-energy £(3,4) so far 
remains in the form given in the Appendix A, Eq. (A5) 
which depends on the auxiliary potentials. In the next 
two Subsections we give self consistent equations for the 
phononic and electronic part in terms of the Green’s func¬ 
tion and the vertex, thereby removing this dependence on 
the auxiliary potentials. 


B. Phonon Propagator, Vertex and Self-Energy 

In this Subsection cast the functional derivative with 
respect to the auxiliary potentials into a dependence on 
the Green’s function only. As a preparatory, we introduce 
the bare phonon vertex 

r° h (ll,22;qT 3 ) = (vi)j - 2 gi(riq)b 12 b T3Tl , ( 17 ) 

where v = ( r~<Jo T z a x tq<j v t z <j z ) and we use the ab¬ 
breviation a q =bq + bL q . The self-energy £ of Eq. (A5) 
in the Appendix A consists of the parts £(1,2) = 
£ xc (l,2) + £f h (l,2). Here, £ xc (l,2) = £ c (l,2) + 
£ ph (l, 2) which separates an electronic and a phononic 
part. Thus, phonon self-energy contributions to £(1,2) 
are given by 

£ H h (l,2) = r° h (l,2-,qT 3 )(a q (T 3 )) T ( 18 ) 

£ ph (l,2) = r° h (l,3-,qT 4 )-G(3,4)- 

■r ph (4,2;q , T , )D ph (qT 4 ,-q , T 5 ) . ( 19 ) 

The Hartree part £ ph (l, 2) is separated in order that 
we obtain a phononic vertex equation with a bare term 
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Eq. (17), later. We have introduced the phonon propaga¬ 
tor D ptl (qr, q'r') = 5 (a- q ' (r')) T / bJ q (r) and the phonon 
vertex r ph (l , 2; qr) = bO~ 1 (l, 2)/b(a q (r)) T . Evaluat¬ 
ing commutators with the Hamiltonian Eq. (5) and using 
fi q = £ 2 -q we obtain the equation of motion 

Td' 2 T a q (r) = 2Q 2 a q (r ) 

- 2n q r° li (ri,22-,-qT 3 )T#(ri)4'\22 ). 

( 20 ) 


Further using Eq. (A4) it is straight forward to derive a 
Dyson equation for D ph (qr, q'r') 


D ph (qr, q'r') = D 0 ph (qr, q'r') + D 0 pii (qr, kr k ) X 
x r° h (ll,22-,-qT 3 )0(33,22) x 
x G(ll,44)r ph (44,33-k'T(,) x 
X D pil (k'T' k ,q'T') , (21) 

with the bare phonon propagator 

(d 2 - f2 2 ) D ph (qr, q'r') = 2f2 q 5 q ,_ ql b(r - t') . (22) 


The phonon vertex in turn satisfies 


r ph (l,2-,q , T l )=r° h (l,2-,qT 3 ) + 


5 S XC (U,22) 


phvi"!yoy ■ 5 0 ( 44 , 55 ) 

xO(44,66)P ph (66, 77; qr 3 )0( 77, 55). 

(23) 


Thus all terms in the Eqs. (18) and (19) do not dependent 
explicitly on the auxiliary potentials. 


C. Coulomb Vertex and Self-Energy 

Repeating the above procedure for the electron- 
electron interaction leads to an expression for the elec¬ 
tronic self-energy in terms of the electron vertex and the 
screened interaction. Using the chain rule, the electronic 
self-energy in Eq. (A5) can be cast into 

S c (l,2)=w(l,4)r z a 0 -G(l,3)-r c (3,2-4). (24) 

We introduce P c (l,2;3) = —SO~ 1 (l,2)/8<P(3) as the 
electronic vertex and derive its Dyson equation with the 
notation P°(l,2-,3) = <Jqt z 8 1 ^ 1,3 

x G(44,66)r c (66, 77; 3)G(77, 55 ). 

(25) 

At this point we have a Dyson equation for the electronic 
single particle Green function (Eq. (14)), the phonon 
propagator (Eq. (21)) and the vertices (Eqs. (23) and 
(25)) to construct the self-energy self-consistently. The 


+ 





+ it 



Figure 1. Diagrammatic representation of the Hedin equa¬ 
tions for a superconductors. 


self-energy E c contains, in addition, the screened interac¬ 
tion w(l , 2) = w 0 (l,3)b$(3)/bip(2). Inserting the def¬ 
inition of the Hartree field and using the chain rule we 
conclude that w(l, 2) satisfies 

w(l,2) = w 0 (l,2) + w 0 (l,3)P(3,4)w(2,4) (26) 

where the polarization P(l,2) = bn(l)/b<P(2) can be 
cast into 

P(l,2) = iw na {a 0 r z • 0(1,3) ■ r c (3,4;2) -0(4,1)}, 

(27) 

where we have used n(l) = Tr na }aoT z ■ 0(1,1)}/2 which 
concludes our self-consistent Hedin cycle. This collection 
of self-consistent equations is represented diagrammati- 
cally in Fig. 1 . In P(l,2), there appear contributions 
from the Cooper channel that add to the screening of the 
bare Coulomb interaction according to Eq. 26. For exam¬ 
ple, inserting the bare Coulomb vertex for P c in Eq. 27 
the term —\f d5Tr a {G(I, —1, 3,1)0(3, 1, 1, —1)} with 
the SC Nambu off diagonal parts of the Green’s function 
contributes to the polarization and thus directly to the 
screening. 

In practice, often gi(rq) is taken as 
’Yhbv eci ( r , h/j!)/ bu q of a phonon calculation 
in density functional perturbation theory (DFPT) and 
similarly the phonon mode frequencies are calculated 
from the density response of a Kohn-Sham system in 
density functional theory 23 . The fl q obtained in this way 
include electronic screening already and often compare 
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well with experiment. In this case, care has to be taken 
not to double count electronic screening, when Eq. (21) 
is applied 24 . 

The question may be reformulated to what phonon sys¬ 
tem with frequencies ft q as “bare” we should start from in 
order that a coupling to the electrons via Eq. (21) leads 
to the best dressed phonon frequencies. While for calcu¬ 
lating SC it is often enough to neglect the influence of the 
electronic system on the phonon system beyond DFPT, 
it would be clearly interesting to investigate the influence 
of the SC phase on the phonon frequencies. This issue 
must be regarded as an open problem as of now and a 
further discussion is beyond the scope of this paper. In 
the following we shall rely on the usual approach to con¬ 
sider the phonons as bare D ph —> D° h and furthermore 
we follow an argument of Migdal 25 and Eliashberg 5 and 
do not consider vertex corrections, i.e. T ph —> P° h . 


IV. THE CONTINUITY EQUATION 

The external pairing field Zl oxt and /l* xt , as well 
as their respective self-energy renormalization on the 
Nambu off diagonal, act as a source and sink for elec¬ 
tron pairs, respectively. Consequently, these terms alter 
the conservation laws obeyed by the exact system as com¬ 
pared to the normal state system. Based on the Dyson 
Eq. ( 14 ) in the form G" 1 ^,;?) = G~\l,2) - £(1,2), 
we generalize the approach of Baym and Kadanoff 26 to 
superconductors. In particular, we derive the continuity 
equation for the electronic charge. First, we realize 

6 1 ^=G(1,3)-G~ 1 (3,2) 

- G(1 ,5) • {2(3, 2) - $(8)<T 0 T z 5 9 ,e) (28) 

5 I , Z = G~ 1 (1,3)-G(3,2) 

- (E(l,3)-$(l)aoT z 5 1 , s ) -G(3,2) (29) 

where (-t 0 <7o3 Ti 6j i3 - H 0 (l, 3))G 0 (3, 2) = f> 1<2 and 
thus 

G-\l,2) = (-T 0 aod Tl 5 1>2 ~H 0 (l,2)) ( 30 ) 

Noting that the derivative is anti-symmetric, 
i. e. d Tl bi t g = -d T2 b lt 2 and V ri 6i,g = -V ri 5 ii2 , we 
calculate 

j d 3(G- 1 (1,3)-G(3,2)-G(1,3)-G- 1 (S,2)) = 

= Jd3(Z(l,3)-${l)a 0 T z 6 1>3 ) -G(3,2) 

— G(1 , 3) ■ (£(3, 2) - $(3)ct 0 t z 5 3 , 2 ) ( 31 ) 

in the Appendix B. Similar to Baym and Kadanoff 26 , 
in addition taking into account the presence of the SC 
condensate many conservation laws can be deduced from 
this equation. We are particularly interested in the local 
2 1 limit of the 1,1 Nambu component of Eq. (Bl). If 


we take the trace in spin space, we arrive at the continuity 
equation for the electric charge in a SC 

i d Tl (h(l)) + V ri ■ ( j(l )) = 


= ^(zl(l,3)-Tr e {&G(3,-l,l. 

. 1 )}) 


- ^[£l(l , 3) ■ Tr B {<frG(3,-1, 1 ,1)}) . 

(32) 

Here, we have used 



(n(l)) = Tr s {a 0 G(l, 1 , 1 , 1 )} 


(33) 

(j-(l)>=Tr> 0 (-(V ri - V r2 ) + 


+ A cxt (l))G(l,l, 2 ,l)| 2 _ 

> 1 } 

(34) 


The right hand side of Eq. (32) requires some further 
explanation. To arrive at Eq. (32) we have used that 
G(3,— 1,1,1) = (G(i,l,S,-l) T 0* an d introduced the 
four component vector (i = 1, 2, 3,4) 

-Arti(l,2) (35) 

Zl i (l,2)=±Tr.{*J'3 m {2(l,l,2,-l)}}. (36) 

In this relation, we use the adjoined of a time and 
space and spin and Nambu matrix as (£(1,2))^ = 

(£(2,l)*) Tlia . Thus, we see that all source terms for 
electronic charge depend explicitly on the SC Nambu 
off-diagonal terms and vanish in the normal state as ex¬ 
pected. The self-energy is not hermitian in general. We 
may, however, decompose it into hermitian and anti her¬ 
mitian parts by 

« t{£(l,2)} =^(£(1,2) +£(1,2)1) (37) 

3m{£(l,2)} = - £(1,2)') (38) 

and thus by definition 

me{£(l,-1,2,1)}= ((fKc{r(2 ,l,l,-l)}) Ts )* (39) 

3m{£(l , -1, 2, 1)} = - ((2fm{r(^, 1, 1 , -1)}) Tb ) (40) 

The four component vector (£'!\) is a spin sym¬ 
metrized composition of the (anti) hermitian self-energy 
and external pair potential. Any spin matrix can be de¬ 
composed into spin symmetric and antisymmetric com¬ 
ponents according to = Tr B {3? Ts A}/2 • 3*^/. This 
is because the Pauli matrices are traceless but cr 2 = cfq 
so that Tr a {(JiA} is a projection on the at axis of a spin 
matrix A. Thus Tr H {$ ls i7}/2 symmetrizes the super¬ 
conducting Nambu off diagonal part of the self-energy 
into spin anti symmetric (4>i) and symmetric ( 4 > 2 , 4 * 3 , 4 > 4 ) 
parts. In a static scenario d Tl (h(l )) = 0 while in the time 
dependent case Eq. (32) remains essentially unchanged 
except that i d/dri —> d/dz where z is on the Keldysh 
contour 22 . 
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V. THE EFFECTIVE INTERACTION 

Among the most important questions in the context 
of unconventional SC is where spin and other electronic 
fluctuations appear in the formalism. The wording fluc¬ 
tuations is often used meaning beyond mean field theory. 
In mean field theory any two particle operator is replaced 
with a single particle operator that couples to the aver¬ 
age in the resulting approximate system of the respective 
other particle. In our case, where the electrons inter¬ 
act via the Coulomb potential and with a phonon field, 
the mean field theory is Hartree-Fock. Similarly, BCS 
theory 2 ' is a mean field theory 18 , however not on the level 
of the Coulomb interaction but on the level of an effective 
interaction that couples single particle states. This effec¬ 
tive interaction may be based on e.g. spin flip processes 
which cannot be described in a mean field theory of para¬ 
magnetic electrons that interaction via a Coulomb poten¬ 
tial. As noted in the introduction Sec. I we distinguish 
fluctuation contributions in our self-energy as the con¬ 
tributions beyond the screened Coulomb diagram. The 
reason is that the unscreened Coulomb potential w 0 does 
not appear in the self-energy and is “very far” from the 
screened potential in a metal. Thus, first, we separate 
the screened Coulomb diagram from the rest of the elec¬ 
tronic self-energy and continue the discussion based on 
the terms beyond the screened Coulomb diagram. We 
generalize the approach of Essenberger et al. 10 and in¬ 
troduce the particle-hole propagator A p for a SC by 

r G (l,2-,3) = r°(l,2-,3) + A p (ll,22,U,55) x 
x G{44,66)r°{66, 77; 3)G{77, 55 ), 

(41) 


where A p satisfies the Dyson equation 


A p (l 1,22,33,44) 


8S{11,22) 8 E{ll,22) 
~ 8 G(33,44) + 8G{55, 66) 
x G{55, 77 )G{88, 66) x 


X 


x A p {77,88,33,44). (42) 


This may be written diagrammatically as given in Fig. 2 

a) . Inserting Eq. (41) into the electronic self-energy 
Eq. (24) our separation of higher vertex corrections in 
terms of the particle-hole propagator is complete. The 
procedure is represented using diagrams as given in Fig. 2 

b) . As an important fact we realize that the screened in¬ 
teraction w{l,2) Eq. (26), similar to the bare Coulomb 
interaction in Hartree-Fock, is independent on spin and 
the bare vertex conserves angular momentum. Thus, the 
first self-energy diagram in Fig. 2 b) cannot account for 
spin flip processes that are discussed as the possible pair¬ 
ing mechanism in the context of High-T c superconductiv¬ 
ity. The screened Coulomb exchange diagram will, how¬ 
ever, contain part of the charge-fluctuations in the sense 
of the fluctuation-exchange approximation 28 . We con¬ 
clude that angular momentum transfer processes are de¬ 
scribed by the higher order vertex contributions, i.e. the 






Figure 2. a) Diagrammatic representation of the Dyson equa¬ 
tion for the SC particle-hole propagator, b) The electronic 
self-energy in terms of the particle-hole propagator, c) Inter¬ 
pretation of the vertex part as an exchange diagram with an 
effective interaction as indicated by the box. d) Contributions 
to the self-energy from the effective interaction I. The arrows 
indicate that the electrons may have different spin. 


second diagram of Fig. 2 b). This second self-energy con¬ 
tribution can be twisted to look like an exchange diagram 
in terms of an effective interaction as indicated by the 
thin box in Fig. 2 c). In the following we refer to the ef¬ 
fective interaction in the box as I {11 , 22, 33 , 44) where 
we have labeled the degrees of freedom 11,22,33,44 
starting with 11 at the left bottom of the box and then 
continuing clockwise. 

This effective interaction has very interesting proper¬ 
ties. For example it can propagate angular momentum 
and other quantities that are conserved at the bare vertex 
due to underlying symmetries. This is possible since the 
diagrams that constitute the effective interaction contain 
a subset that has no Green’s function connection between 
“top” {22,33) and “bottom” {11,44). 

An example of such contributions are the horizontal 
ladder diagrams beyond second order. First and second 
order are part of the screened Coulomb self-energy and 
thus have to be excluded to avoid double counting. In 
this type of diagrams, the conservation at the bare vertex 
appears independently on the connected top {22 t— 33) 
and bottom {11 —> 44) Green’s function lines. Instead, 
top and bottom are connected only with interaction lines 
and thus the diagram, viewed as a whole, does not obey 
the conservation constrains of the bare vertex. Examples 


















7 


of such contributions are given in Fig. 2 d). For the 
second contribution we have chosen the diagrammatic 
arrangement of Doniach and Engelsberg 29 . 

For example, a paramagnetic system subject to a 
Coulomb interaction has a bare vertex that is propor¬ 
tional to (Jo in spin space. Still, this type of effective 
interaction in the box of Fig. 2 c) can propagate an¬ 
gular momentum in the sense that it has components 
proportional to <J x ,y,z as well, i.e. it describes spin flip 
processes. Similarly, it also allows to describe how off di¬ 
agonal Nambu components are rotated due to the inter¬ 
action. The bare electron-electron and electron-phonon 
vertex in Nambu space is proportional to t z instead of a o 
as for the spin space, due to the fermionic commutation 
rules. The box of Fig. 2 c) corresponds to an interaction 
that has also components in t X)V and thus in the language 
of Anderson 1 ' propagates a rotation of iso-spin. These 
components correspond to the Cooper channel and an in¬ 
teraction mediated by the fluctuation propagator in the 
sense of Ref. 14. 

1(11,22, 33,44) is a function of four N ambu and spin 
indices and time and space variables, respectively. We 
suppress the time and space arguments and decompose 
1(1,2,3,4) along the Nambu and spin basis 

I(1,2, 3, 4) ='52 I ° J j((TaTi®<J b T j )- 1 2S4 (43) 

ijab 

!ij = I C 1 ^JA)(^aTi)i2((JbT J ) n (U) 

1234 

We write (t* ® Tj)j 234 = ( T i)l 2 ( T j )34 to indicate that 
the basis vectors t, and r ; - form an outer product. In the 
following we discuss only the Nambu degrees of freedom 
since the spin is completely analogue. Considering the 
Nambu x,y symmetric part of the effective interaction, 
we may further decompose 

y, IijTi ® Tj = I-+(t x - i Ty) <g> (t x + i Ty) + 

ij=x,y 

+ I+-(r x + i Ty) <g> (r x - i Ty) + 

+ I++(r X + i Ty) (g) (r x + i Ty) + 


with 

+ / - (T x - i Ty) <S> (T x - i Ty)) 

(45) 


I — b — {Ixx “1“ lyy H-xy H - i lyx) !4 

(46) 


I -\— — {Ixx “1“ lyy 4“ i Ixy i lyx) !4 

(47) 


-^++ = (Jxx lyy 4“ H-xy 4“ i lyx)/ 4 

(48) 


I -— [Ixx lyy H-xy Hyx') / 4 

(49) 

Note that for any complex number Z we find 


(t x 

± iT y )® = sign(3?®)|®|(r a; ± ir y )e TIT ^ 

(50) 

(■ T v 

± ir 2 )® = sign(3?®)|®|(r y ± i T z )e ¥1T ^ 

(51) 

(Tz 

± i t x )Z = sign(3?®)|®|(r z ± i T x )e T1Ty,p 

(52) 


where <f> = arctan(3.21/3?®). Also note thatT r e =FlTj! ' ! '^ = 
e ±iT z . y 4> Tx anc j s j m }i ar f or other Pauli matrices. Thus, 


we may write 

(t x ± i T y ) ® (t x ± i T y )I±± = sign(3?/±±)|/ ±± | x 

x ( t x ± ir y )e =FIT!, ^ ± <8 
® e ±1T * (t x ± i T y ) (53) 

where 

<f>±± = arctan(3/±±/3?/±±). (54) 

Since the effective interaction I is used in an exchange 
diagram, we note that (t x ±ir y )•G• (r s ±ir y ) interchanges 
the (±1, =Fl) component of the Green function with the 
(=Fl, ±1) component in Nambu space. The (±1, =Fl) com¬ 
ponent is the complex conjugate of the (=pl,±l) com¬ 
ponent. Such terms I±±, constructed from /“ y ,/ y *,/ yy 
and /“£, are the SC analog of a spin flip interaction 
such as the one of Essenberger et al. 10 . Furthermore, 

. <f>±± - _i_. <f>±± 

* ■G-e ±1T * corresponds to a phase rotation of 
the Nambu off diagonal parts by ±<j>±± . Thus, this part 
of the effective interaction corresponds to a (non-local) 
phase transformation of the condensate where the phase 
is rotated by ±</>±± and then conjugated. 

In our notation, the spin flip interaction corresponds 
to the terms I z l J , ///, /// and /// since they correspond 
similarly to a self-energy contribution ~ (a x ± icr y ) • G ■ 
(a x ± i a y ). 

VI. SUMMARY 

We have derived Hedin’s equations for superconduct¬ 
ing electrons that are coupled to a system of phonons. 
Our coupled equations are formally exact. Following the 
approach of Baym and Kadanoff, we derive the continuity 
equation for a SC where the external pair potential and 
its self-energy renormalization appear as source and sink 
terms for electronic charge. We point out how we can de¬ 
fine an effective interaction that can describe fluctuations 
beyond the screened Coulomb diagram, e.g. of spin of the 
electrons and phase of condensate. This interaction does 
not suffer from double counting problems and is given 
rigorously in terms of the particle-hole propagator. 
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Appendix A: The Electronic Dyson Equation 

In this Appendix, we complete the derivation of the 
Dyson equation Eq. (14). We collect the contributions to 
the Dyson equation of the single particle parts H 0 + H aux 
of the total Hamiltonian Eq. (5) on the left hand side 
of the equation of motion Eq. (13). The commuta¬ 
tors T[IJ 0 _ e + H._ V ,^}(1 y)\ and T[7F e _ e + (1 y)\ 


are straight forward to evaluate component-wise while 
the time ordering symbol is essential to cast the result 
back into the unified notation with the Nambu operators 
&(11). We obtain 


E k3(s ljS (TQ(7od Tl + a 0 T z tp(l))H 0 (l, 3)^ __ x 

3 J 

x G(33, 22) = 

= ^E Jd3w 0 (l,3)(aoT z ) S} s(aoT z )jj x 
x G w (33,11,33,22) 

- E (vi)- 1 3g i (r 1 q)(a q {T 1 )<F{13)^(22)) T . (Al) 

3,q,i 


with the two particle Nambu Green’s function 

G (2) (l 1,22, 33,44) = {&{! 1)&(22)$(S3)^(44)) T . 

Furthermore we derive the relation 


TT aa {- 2 a 0 T z -G(3,3)}G(ll,22) 


1 


Furthermore, for a generic field A 

(A4) 

At this point we introduce the self energy 

2 ( 1 , 2 ) = 


Wq( 1 , 3 )t z <jq -G(l,4) 
-r° h (l,3;qT 4 )-G(3,4) 


bG-\4,2) 

M3) 

6G~ 1 (4,2) 

s Jq{r 4 ) 


+ 


+r°(l , 3\qr 4 ) ■ 6 2i3 croro(M(T 4 )). 


(AS) 


= E^Mo r z y 3t - s G^(33,11,33,22) 


(A2) 


Inserting the Eqs. (A2) and (A3) into Eq. (Al), together 
with the definition of the Hartree Hamiltonian H h (1,3), 
we arrive at 

Jd3(-6 1 , 3 T 0 a 0 d Tl -H H (l,3))-G(3,2) 

= 6is6n+JdSE(l,3)-G(S,2). (A6) 

Now we insert GhG” 1 and with the Hartree Green’s func¬ 
tion Eq. (15), applying G H from the left, we arrive at 
Eq. (14). 

Appendix B: Conservation conditions for the Green 
Function 


and similarly 


8G(11,22) 


6J 9 (r) 

<o,(ti )£(1 ~ 1 )&( 22 )) t . 


— (a q (Ti)) T G(l 1,22) = 


(A3) 


In this Appendix, we want to discuss the conditions 
implied by the Dyson equation for the electronic Green 
function Eq. (14). We give the non-local basic equa¬ 
tion, that we use to derive the continuity equation. The 
procedure is a straight forward adaption of Baym and 
Kadanoff 26 , noting that G” 1 and G are 4x4 matrices 
that do not commute. The result for the left hand side 
of Eq. (31) is 


J d 3(G-\l,3).G(3,2)-G(l,3).G-\3,2)) = 

= (^°o T o(d Tl + d T2 ) — o'ot z -(— i( V ri +V r2 ) + A ext (l) — A oxt (2)) (—i( V ri — V r2 ) + A oxt (l ) + A oxt (2 )) 
+ c 0 t 2 (</> ext (l) - t(«)) + * (S(r 0 + t z ) - S*(to - Tz )) ■ (B^l) - £ cxt (M)) -G(l,2) 

- ((iV ra + A cxt (2)) 2 + </> ext (2)) [a 0 r z ,G(l , 2)] - ^[(S(r 0 + r z ) - S*(r 0 - r z )) ■ B ext (2),G(l, 2)] 

+ \j dS((T x -iT y )AUl,S)-* - (T x +ir v )At«(l,S)-*) - G(3,2) 

~\j d 3G(1,3)- ((t x -\.t v )A^( 3,2)-*-(T x + iT y )A:^(3,2)-*) 


(Bl) 
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The Nambu off diagonal right hand side of Eq. (31) is then evaluated to be 

\ I dS((T x + iT y )S(l,l,S,-l) + (T x -iT y )S(l,-l,S,l))-G(S,l) 

-\j d 3G(1,3) ■ ((r x + i Ty)S{S, 1, 2, -1) + (r x - i t v )S(S, -1, jS, 1)) 

_ f 3m{/ dSfRt{H(l,l,S,-l)}G(S,-l,2,l)} 91 e{f dSXe{ij(J,l,S,-l)}G(S,-l,2,-l)} ) 

~ V 9R{/d3iRe{r(I,-l,3,l)}G(3,l,£,l)} 3m{fdSiHe{£(l,-l,S,l)}G(S,l,2,-l)} J 

( 9te{/ d33m{S(l , 1, 3, -1)}G(3 , -1, 2, 1)} 3m{/ d33m{E(l , 1, 3, -1)}G(3, -1, 2, -1)} \ , , 

+ ^ 3m{/d33m{r(l,-l,3,l)}G(3,l,S,l)} JKe{/ d33m{£(l ,-1, 3,1)}G(S, 1, 2,-1)} / ' 1 j 

The Nambu diagonal part can be computed in an analogous way. Taking the local limit 2 —> 1 of the Eq. (Bl) gives 


J d 3(G- 1 2 3 4 5 6 7 * 9 10 (1,3).G(3,1)~G(1,3).G~\3,1)) = 

= 2itoto<9 T i G(1, 2)\^->-i + 2<roT z iV T . 1 ( — (V ri — V r2 ) + A cxt (l ))G(1 , 2)\%^i 
- ((iV r2 + A oxt (l )) 2 + (j) ext (l [a 0 r Z: G{l , 2)]\g->i - ^[(S(r 0 + t z ) - S*(r 0 - r z )) ■ G(1,1)] 

+ IJ dS((r x +iT v )A^{l,S) ■ $ - (t x - W v )Al t (l,3) ■ d?) • G(3,l ) 

~IJ d3G(l,8) - ({T x + iT v )A^{S,l)-*-(T x -iT y )AZ«(S,l)-*) . (B3) 


It is important that the commutator [uot z ,G(3, 2)] has 
only components on the Nambu off diagonal. If we are 
interested in the usual conservation laws of electronic 
charge or magnetic density, these terms do not appear. 
Taking the local limit and selecting only the 1,1 compo¬ 
nent, the equations simplify significantly and we arrive 
at Eq. (32). 


Furthermore, we point out that if the self-energy 
is non-SC, i.e. has only Nambu diagonal non van¬ 
ishing componentsi7(l, 3) oc t z ,t 0 , the right hand 
side of Eq. (31) is proportional to the commutator 
[<tot~, G(3, 2)\. This commutator on the other hand does 
not have components on the Nambu diagonal, i.e. the 
usual continuity equations for electronic charge and mag¬ 
netic density are satisfied without self-energy contribu¬ 
tions as expected. 
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